Dynamic Statistical Scaling in the Landau-de 
Gennes Theory of Nematic Liquid Crystals 



By Eduard Kirr, Mark Wilkinson and Arghir Zarnescu 



Abstract 

In this article, we investigate the asymptotic behaviour of a correlation function 
associated with a nematic liquid crystal system undergoing an isotropic-nematic phase 
transition following an instantaneous change of system temperature. Within the 
setting of Landau-de Gennes theory, we confirm a hypothesis in the condensed matter 
physics literature on the average self-similar behaviour of the correlation function in 
the asymptotic regime at time infinity, namely 

H 2 _i 

"u (r,t)-e 8t =0(t 2) as t — > oo, 

L°°(R 3 , dr) 

where the correlation function c M0 : R 3 x (0, oo) — > R is defined by (6) below. In 
the final sections, we also pass comment on other possible scaling behaviour of the 
correlation function. 



1. Introduction 

Nematic liquid crystals form a class of condensed matter systems whose constituent 
rod-like molecules give rise to rich macroscopic nonlinear phenomena. Thermotropic 
nematic liquid crystals are a well-studied subclass of nematics whose physical properties 
change dramatically with variation of system temperature. It is commonly observed 
in the laboratory that above a certain temperature threshold, the rod-like molecules 
exhibit no local orientation preference, which is known widely as the isotropic phase of 
the material. Yet, reducing the temperature of the material below this threshold results 
in the molecules arranging themselves along locally preferred directions, yielding the 
so-called nematic phase. 

Models for the dynamics of this transition from isotropy to a nematic phase in 
liquid crystals present many challenging mathematical questions, both in continuum 
theories and also in statistical mechanics. In the laboratory, the transition between the 
isotropic and nematic phases is observed to be effected by the seeding and subsequent 
growth of nematic 'islands' in the ambient isotropic phase. The characteristic length 
scale L = L{t) of these structures increases with time once the system is coerced into a 
transition of phase. It is argued in the physics literature that the domain coarsening of 
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the nematic phase is a scaling phenomenon, namely the structure of the nematic profile 
at late times looks statistically self-similar to those profiles at an earlier time, up to a 
simple rescaling in space. We invite the reader to consult the extensive review article of 
Bray [4] for a helpful introduction to this area of condensed matter physics. Laboratory 
experiments performed by Pargellis et al. [16], in which a nematic liquid crystal 
system was contrived to resemble a dynamic XY- model, support this argument. 

It is the aim of this paper to articulate such behaviour of thermotropic nematic 
systems in a mathematically rigorous manner. We work within the framework of con- 
tinuum PDE theory, as opposed to the forbidding setting of statistical mechanics. 

As the isotropic-to-nematic phase transition is inherently dynamic, we require an 
appropriate evolution equation to model this coarsening phenomenon. The model we 
employ in this paper is the L 2 -gradient flow of the well-studied Landau-de Gennes 
energy. This energy is given by 

EhdG [Q] ■= / r3 Q|VQ| 2 + ^tr [ Q 2 ] - |tr [Q 3 ] + \ (tr [ Q 2 ])^ dx (1) 

for appropriately regular maps Q, the gradient flow of which is 

^ = AQ - a{9) Q + b(Q 2 - itr [Q 2 ] /) - ctr [Q 2 ] Q. (2) 

We invite the reader to consult the works of Ball and Zarnescu [2], Ball and 
Majumdar [3], Majumdar and Zarnescu [13] and Nguyen and Zarnescu [15] 
for an introduction to analytical and topological aspects of Landau-de Gennes theory. 

In order to achieve our goal, we must first state in precise quantitative terms what 
we mean by structure of nematic profiles and also by statistical self-similarity. Indeed, 
both of these concepts lead us naturally to the mathematical problem of obtaining a 
phase portrait for a quantity (a correlation function) describing the average behaviour 
of solutions of the gradient flow of (1) above. We argue that one might view this 
as a convenient simplification of the infinite-dimensional dynamics generated by this 
PDE on its natural phase space, which captures the essential qualitative behaviour of 
solutions. However, antecedent to both of these ideas is the notion of order parameter, 
which we now discuss in section 1.1 before introducing correlation functions in section 
1.2 (related to structure) and statistical solutions of PDEs in section 1.3 (related to 
statistical self-similarity). 

1.1. The de Gennes Q-tensor Order Parameter. The first task when building 
a continuum model of such phase transition phenomena is to decide upon an appro- 
priate order parameter that captures small scale material structure and allows one to 
distinguish between different phases of the material under study. In this article, we 
work with the Q-tensor order parameter which is able to describe both uniaxial and 
biaxial phases of nematic liquid crystals, as opposed to, say, the director-field formalism 
of Ericksen-Leslie theory [11] or Oseen- Frank theory [9]. 

Let us suppose that to each point x in a material domain Q C M 3 we associate a 
probability density function p x on molecular orientations which lie in § 2 . In order to 
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model the Z2 head-to-tail symmetry of nematic molecules, each density is endowed with 
the antipodal symmetry p x {n) = p x (—n) for all n £ S 2 . It was the idea of the Nobel 
prize-winning physicist Pierre-Gilles de Gennes that one might consider the essential 
macroscopic information of the system to be contained in the matrix of second moments 
of p x with respect to molecular orientations. The de Gennes Q-tensor order parameter 
(with respect to the density p x ) is defined to be 



It is a normalised matrix of second moments of p x , and one may quickly check that 
Q(x) is a traceless and symmetric 3x3 matrix. 

At the level of probability measures on molecular orientations, the uniform density 
p = 1/|§> 2 | corresponds to the isotropic phase of nematics. It is important to note that 
the term —1/3/ (which contains no information about the system) is included in the 
above definition (3) of the Q-tensor by convention, so as to render Q identically zero 
in the isotropic phase. The Q-tensor may then be thought of as a macroscopic order 
parameter that measures the deviation of the system from isotropy. For a more detailed 
introduction to the Q-tensor order parameter, one might wish to consult de Gennes 
and Prost [5], Majumdar [12] or Newton and Mottram [14]. 

We denote the order parameter manifold of all such matrices by Sym (3), namely 



where tr : M 3x3 — y R denotes the familar matrix trace operator. We subsequently refer 
to maps which take their values in Sym (3) as nematic profiles. 

1.2. 'Structure' of Nematic Profiles: Correlation Functions. The measure 
of structure of nematic profiles Q : K 3 x (0,oo) — > Sym (3) we employ in this paper is 



One may think of c as quantifying how correlated different regions of a nematic profile 
are with one another and thus as some measure of the spatial structure of the system 
as it evolves over time. For further information on the use of correlation functions 
in condensed matter physics, we recommend that the reader consult Sethna ([19], 
Chapter 10). 

As mentioned above, it is argued in the physics literature (see for instance Den- 
niston et AL. [6] or Zapotocky et al. [21]) that once a system is coerced into an 
isotropic-nematic phase transition the correlation function assumes a self-similar form 
asymptotically in time. In summary, 




(3) 



Sym (3) := {Q G M 3x3 : Q T = Q and tr[Q] = o} 



the two-point spatial correlation function c : M 3 x (0, oo) — > M defined by 




(c(r,t)} ~r 



( 



L(t) 



r 



) 



in an appropriate topology as t 



oo 



(4) 
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where T is some 'universal' scaling function and L(t) is the aforementioned character- 
istic length scale of nematic domains which invade the isotropic phase. The angular 
brackets (•) signify that one considers a suitable average value of the correlation func- 
tion as computed over many repeated experimental trials, smoothing out anomalous 
data gathered from experiment. If one then wishes information on the structure of this 
system for large values of time t, one need only rescale in the spatial variable r. 

In order to model such an averaging procedure in a mathematically rigorous man- 
ner, we employ the notion of statistical solution of evolution equations introduced by 
FoiA§ [7] and also by Vishik and Fursikov [10]. This notion of solution to nonlin- 
ear PDE was subsequently studied many authors in the mathematical theory of the 
Navier-Stokes equations to explore aspects of turbulence theory. We now outline the 
significance of statistical solutions in our study of phase transitions. 

1.3. 'Statistical Self-similarity': Statistical Solutions of PDEs. One can 

show that the semiflow {S(t)} t >o of the gradient flow (2) of E^g possesses neither 
a global attractor nor an inertial manifold in the natural phase space 
L 6 (M 3 ), which is unfortunate as both objects are well suited to revealing asymptotic 
structure of dynamics. Since we are interested in studying the time asymptotics of the 
correlation function c^ , this presents a significant challenge. Thus, acknowledging the 
structure of the manifold of solution trajectories 



Qo&H <- > 
as being rather complicated, we ask naively what one could say if one turns instead to 
the manifold of measure-valued statistical solutions 



where Mq(H) denotes the set of all Borel probability measures on H. The manifold 
M is 'coarser' than A4 in terms of information, in the sense that each such statistical 
solution nt of (2) at time t is insensitive to changes on sets of measure 0. In order to 
make use of this insensitivity to 'fine structure' of the trajectory manifold A4, and to 
perhaps also give us hope of characterising one qualitative aspect of solutions of (2) on 
H (section 6), we redefine the correlation function c as 



It is with this definition of the correlation function we work throughout this paper, as 
opposed to the earlier- given definition (4) which is supported only on single trajectories 
in Ai. Moreover, it allows us to give good mathematical sense to the average (•) 
discussed above, namely (c(r,t)) := c^ (r,t). 

The gradient flow (2) takes the shape of a nonlinear heat equation for the evolution 
of nematic profiles Q. As this equation generates a strongly-continuous semigroup of 







(6) 
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solution operators {S(t)}t>o on phase space H , it is possible to deal more concretely with 
the evolution of initial measures as opposed to the 'Liouville equation' formalism 1 of 
statistical solutions, which can be adopted when a semigroup of solution operators is not 
readily available. To be precise, for a given suitable initial Borel measure /io supported 
on H, we construct in Section 4 an associated one-parameter family of measures {nt}t>o 
with the property that 

= ^ E ^Awq,, ( E ) 

i=i 

for measurable subsets E C H, nematic profiles Qj t k £ H and t > 0. One can then 
see each statistical solution {/J.t}t>o £ as a convenient concatenation of individual 
solutions {Q{-,t)}t>o S -M. 

We store some more remarks on the construction of a phase portrait for the corre- 
lation function c^ in the closing section 6. 

1.4. The Mathematical Model. As we work with tensors, we consistently em- 
ploy the Einstein summation convention over repeated lower indices. We now present 
a brief derivation of the models (1) and (2) introduced above. 

In accordance with the Landau theory of phase transitions, the free energy density 
W is assumed to be a smooth map which is frame indifferent, in our case invariant 
under the action of the rotation group SO (3): 

W(T R D,RQR T ) =W(D,Q) 

where Q £ Sym (3) and (T R D) ijk := Y$, m ,n=l RuRjmRknD imn for all R G SO(3) 
whenever D is a rank 3 tensor, i.e. W is hemitropic in its gradient argument and 
isotropic in its field argument. We assume that the corresponding free energy has the 
form 

/ W(VQ(x),Q(x))dx= f (E(VQ(x))+F(Q(x)))dx 

JR3 J R 3 ^ I 

in which the elastic and bulk terms are decoupled. One may show by invariant theory 
that any such smooth isotropic map F is a smooth function Fq : R x R — > R of the 
principle invariants of SO (3) 

tr [Q 2 ] and tr [Q 3 ] . 

We consider the energy functional 

EhdG [Q] ■= [ (\\VQ(x)\ 2 + F (Q(x))) dx, 

to which the elastic energy contribution is simply the Dirichlet energy, although one 
could work with more a general elastic energy of the form 

Li 2 ^ L>2dQikdQij ^ L-jdQjjdQjk ^ ^Q ik ^Qjj 

2 2 dxj dxk 2 dxj dxk 3 dxk dx£ 



See FoiA§ et AL. [8], Chapter V. 



6 



AQ - a{9) Q + b(Q 2 - hsx [Q 2 ] /) - ctr [Q 2 ] Q. (8) 



for Li > 0, which is compatible with the system symmetry. We employ the simplest bulk 
energy density Fq that is able to predict a first-order isotropic-nematic phase transition, 
namely 

W):=^tr[Q 2 ]-^tr[Q3] + |( tr [Q 2 ]) 2 , (7) 

where a = a${9 — 9 C ), 9 C > is the critical temperature and 9 G R is a temperature 
parameter controlling the depths of the nematic energy wells. The material-dependent 
constants (a, b, c) belong to the bistable region 

V := {(a,b,c) G (0,oo) 3 : b 2 > 27ac} . 

of parameter space (0,oo) 3 . The model for the dynamics of nematic liquid crystals we 
study is the L 2 -gradient flow of the functional L^q, 
dQ 

dt _ " " ' " 3' 

We consistently refer to (8) as the Q-tensor equation. Furthermore, we model the 
temperature quench as instantaneous, i.e. the parameters (a, b, c) lie in the bistable 
region T> and do not depend on an evolving temperature field 9. Before the quench is 
effected at time t = 0, the parameters (a, b, c) lie in the complementary region (0, oo) 3 \ 
D: in this case, the isotropic phase Q = is a global minimiser of the bulk energy 
density (7) above. 

We work over the whole space M 3 so we need not worry about restricting the 
spatial domain of the correlation function (6). One may interpret this physically as 
modelling behaviour in the bulk of the nematic material, far from any boundaries where 
complicated boundary effects could come into play. 

The main result of this paper is captured by the following statement. 

Statement of Main Result. For any given 5 > there exist rj > 0, depending 
only on 5 and the parameters (a,b,c) G T>, and an open dense subset 2 of 

R G L°°(M 3 ,Sym (3)) : esssup (1 + \x\) 8+8 \R(x)\ <n 

such that for any Borel probability measure supported in this open dense set, the asso- 
ciated correlation function c^ (6) exhibits asymptotic self-similar behaviour, viz 



_lrf 



Olt 2] as t — > 00. 



L°°(R 3 ,dr) 

1.5. Structure of the Paper. In Section 2, we record some preliminary results 
which allow us to establish a decomposition formula for 'small initial data' solutions of 
the Q-tensor equation in Section 3. The main result of the paper lies in Section 4, in 
which we establish asymptotic self-similarity of the correlation function c^ by means of 
the previously-established decomposition formula. In the final section of the paper, we 
briefly discuss scaling behaviour of the correlation function corresponding to a class of 



The structure of this open dense set is described in section 4.1. 
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solutions starting from 'large' initial data, which gives rise to a different scaling regime 
for the correlation function, namely L(t) = t. 

1.6. Notation. For any matrix A 6 Sym (3) we denote its Frobenius norm by 
\A\ := V A : A, where A : B = AijBji for A,B£ Sym (3). In what follows, unless oth- 
erwise stated, all Lebesgue spaces have range in Sym (3) and we write L P (M 3 , Sym (3)) 
simply as L P (M 3 ) for 1 < p < oo. We also work with the phase space H := L 2 (IR 3 ) n 
L 6 (R 3 ) endowed with the norm ||Q|| := max{||(5||2, ||Q||6}; where the norms || • || p are 
given by 



J (tr [Q(x) 2 ]) 2 dx^ P when 1 < p < oo, 
ess sup |Q(x)| when p = oo. 



We denote the fundamental solution of the heat equation on R 3 by <£>, namely 

e" 



\x\ 2 /4t 



*0M) ■-- rir, 

(4vrt)2 

and denote the time-shifted profile ^(x, t + 1) simply by t). Finally, for any given 
5 > we denote by A the space of essentially bounded functions given by 

A := \q e L°°(1R 3 ) : esssup {I + \x\f +S \Q{x)\ < oo\ , 

and equip it with the natural norm ||Q||„4 := esssup a , gK 3(l + |x|) 8+5 |(5(2 ; )|- The constant 
5 > is chosen arbitrarily but remains fixed. 

2. Some Preliminary Results 

We now establish some results which enable us to comment on the long-time be- 
haviour of the correlation function concentrated on individual solutions of the Q-tensor 
equation, namely the asymptotic behaviour as t — > oo of the quantity 



c(r,t) 



tr [Q(x,t)Q(x + r,t)} dx 



tr [Q(x,t) 2 ] dx 



where the initial datum Qq G H is chosen appropriately. We begin by stating a simple 
result on the well-posedness of the Q-tensor equation in the phase space H. 

Proposition 2.1. For each Qq £ H there exists a unique global classically- smooth 
solution Q of the Q-tensor equation (8) on M 3 x (0,oo) satisfying ||Q(-,t)|| < oo for all 
t > 0. Furthermore, the nonlinear semigroup of solution operators S(t) : H — >• H 
defined by 

{ Qo, if t = 

S(t)Q ■= I 

{ Q(;t;Qo) if t >0 
for any Qq G H is strongly continuous. 
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Proof. The nonlinear semigroup {S(t)}t>o is constructed by standard techniques. 
See Pazy [17] for details. □ 

It is our aim to investigate finer properties of 'small initial data' solutions of (8), 
namely we wish to demonstrate a decomposition result for solutions starting from initial 
data which are small in the || • ||_4-norm. In order to achieve this, we first must examine 
the long-time behaviour of solutions of the heat equation in L P (IR 3 ). 

2.1. A Remark on the Heat Equation in M. d . Consider the heat equation 

^(x, t) = Au(x, t) (x, t) €R d x (0, oo), 

where 1 < p < oo, d > 1 is an integer. For given L q (M. d ) initial data, it possesses a 
classically smooth solution on R d x (0, oo) and this solution u(x, t) may be expressed 
as a convolution involving the d-dimensional heat kernel <&d( x ,t), 

r e -\x-y\ 2 /u 

u(x, t) = ($d(-,t) * u )(x) = / , uo(y) dy, 



Rd (4irt) d / 2 

for t > 0. In order to obtain an estimate on the time decay of solutions in L p (M. d ), one 
can employ Young's inequality in the above to produce 

M / Ml W „ „ 

\\u{;t)\\ p < Tx \\Uo\\q, 

for all t > and 1 + 1/p = 1/q + 1/r. Thus, all solutions of the heat equation with 
L q (M d ) initial data decay to zero at least with rate t~ d / 2+d l 2r as t — > oo. We show in 
what is to come that if we insist that the initial data of the heat equation are of zero 
mean i.e. J Rd uq = and possess some mild decay condition at infinity, then we may 
improve our estimate on the rate of decay of the corresponding solution to in L p (R d ) 
as t — > oo. 

We begin with an auxiliary pointwise estimate on solutions of the heat equation 
that will be of help to us in the following section, and may indeed be of independent 
interest. We consider only the case of the heat equation in M 3 appropriate to our 
considerations. The methods used here may be extended without great effort to capture 
a similar result in M. d for arbitrary dimension d > 1. 

PROPOSITION 2.2. For u e A, consider the difference 

m(x,t) := (e tA u ){x) - $(x,t) / u (y)dy. 

Jr 3 

For any ft > 0, the estimate 

< (i + J|)"7 kj M (i + J|) VoWl *, (9) 

holds for x E M 3 , t > 0. The explicitly- computable constant C > depends only on j3. 



Proof. Define the map a : [0, 2ir) x [0, ir) — > § 2 , a standard surface patch parame- 
terisation on the unit sphere, to be 

<r($, <p) := (cos i? sin <p, sin $ sin <p, cos <p) , 

and define U : [0, oo) x [0, 2vr) x [0, vr) -»• R by 

C(r, i?, <p) := / uo(/?cos #sin<p, psini? sin<p, pcos (p)p 2 sin<^ dp. 

Now, by a change of co-ordinates and an application of integration by parts, we deduce 
that 



(e tA u )(x) 

n r 2n f-oo p -\x-ra(tf,ip)\ 2 /At 



Mo(rcr(^, <p))r 2 sin ip drdfidp 



o Jo Jo (4vrt) 3 / 2 

f d 
= / $(x - ra,t)— (— [/(r-,i?,<p)) drdtfdip 
Jrz or 

= $ T U(0,$,cp)d$dip+ [ ° ' ^ ~ r<7 - ) $(s - r<7, t)I7(r, fl, y>) drd&dtp 
Jo Jo Jr 3 2t 

f f (7 ' [X — T(7) 

= <& / u (y)dy + / -<J>(x - ra, t)U(r, p) drdtfdtp. 

it 3 it 3 2i 

Applying the elementary inequality y / iexp(— z) < exp(— z/2) for z > 0, we deduce 
from the above 

(e tA u )(x) -$(x,t) / u (y)cfy 

* (4^4 Xa e " k_rCT|2/ " ( MHI/^p) sin V drcftW* 
< Cr 2 ( sup ( e -l^-H 2 /s^ [ / >OO | U0 ( pCT )|p 2 dp] smpdrdtidip 

= Ct- 2 J^e^ x ^ r)2/8t (^J^ \uo(pa)\p 2 dp ] j smipdrdtidp, 

where C > is independent of t > 0. A couple of applications of Fubini's theorem then 
yield 

m(x,t)<Ct- 2 f° e -(N^) 2 /8t f \ Uo ( y )\dydr 
Jo J\y\>r 

= Ct ~ 2 L (j^ l -XB{o,r){y))e-^- r)2/St dr) \u (y)\dy, 

where XB(o,r) denotes the characteristic function of the open ball 5(0, r) C M 3 . From 
this we find 

m(x,t)<Ct- 2 [ ( [ M e-fcW** dr) \u {y)\dy. (10) 
Jrs \Jo J 

We pause at this point to obtain a helpful estimate on the y-integrand in the above. 
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Lemma 2.3. For any j3 > 0, there exists a constant C = C(/3) > such that 

J _ e-i di < CY [y^x) for X > Y >°- (11) 

Proof of Lemma. In what follows we fix Y > 0, considering it as a parameter. We 
separate the demonstration of the inequality into two cases. 

(i) The case where < X < Y: One may easily deduce from the chain < X < Y the 
simple inequality 

Y <Y 

for any f3 > 0. Using this, we have 



1 + Y\? 



1 + XJ ' 



X r X f i _(_ y s & 

di < / max (e~^) d£ = Y < Y 1 ' ' 



X-Y JX-Y So ^ 1 + X 

Thus, the claimed inequality holds for < X < Y. 

(ii) The case where Y < X: This time, we have the estimate 

[ X e~?dt<[ X sup(e-«o) d£ = Ye-( x - Y ^. 

JX-Y JX-Y £ 

Now, since there exists a constant C > 1 dependent only on j3 > such that 

1 \? 



e- z2 < C 



l + Z 



for all Z > 0, and noting that 



1 V ( l+Y V 



-(X-Y)* < c 



1 + ZJ ~ \l + Y + ZJ ' 
for all Y > 0, it follows from the substitution Z \— > X — Y that 

i+y_y 

1 + XJ ' 

for all X > Y. This estimate together with the first case estimate yields the result (11) 
above. □ 

We now return to the estimate (10) and close the proof of the proposition. Trans- 
forming the y-integrand through a simple change of variables provides 

fn(x,t)<Ct- 2 jT g (j^e-^-^^dr^ \u (y)\dy 

= Cr 2 [ (V&t P 7 , e-?di ) \u (y)\dy. (12) 

Finally, setting X := \x\/{St) 1 / 2 and Y := |y|/(8i) 1//2 , we apply the above lemma to 
(12), from which estimate (9) follows. □ 

Remark 2.4. Raising both sides to the power p and integrating over the domain 
M 3 , we may use the result in (9) to show that 

<^^/ K(y)|(l + \y\) p+1 dy, 



e tA u -$(-,t) / u (y)dy 
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for t > 0. Therefore, if the initial data of the heat equation are of mean zero and lie in 
L^IR 3 ; duip), where dujp := (l + |x|)^ +1 <ix, we have improved upon the standard estimate 
obtained by Young's inequality for the time decay of solutions in L P (IR 3 ). Estimates 
of type (9) are key to obtaining a decomposition formula in the following section for 
solutions of (8). Let us also remark that the above result leads to a generalisation to 
higher spatial dimension of the work of Taskinen [20]. 

3. A Decomposition Formula 

We establish here a helpful structural formula for solutions of the Q-tensor equa- 
tion starting from 'small' initial data. In particular, we show that the leading order 
asymptotics (up to order i~ 3 / 2 in time after rescaling) of such solutions in all L p -norms 
are governed by a matrix multiple of the heat kernel. It proves important in our sub- 
sequent analysis to consider a transformed version of the Q-tensor equation (8). If we 
set R(x,t) := e at Q(x,t) and substitute R into (8) in a formal manner, we obtain the 
equation 

dtR = AR + b e- at (r 2 - itr [R 2 ] ) - c e - 2at tr [R 2 ] R, (13) 

which is the principle equation of analysis in the section to come. 
Let Xq denote the set of maps 



X := |w € C(R 3 x [0,oo); Sym (3)) 
where the space-time weight u is defined to be 



sup oj(x,t)\W(x,t)\ < CO 



4+f 



When endowed with the natural norm 

\\W\\ Xo ■= sup to(x,t)\W(x,t)\, 

(x,t) 

one may verify that Xo admits the structure of a Banach space. We also employ the 
product space X := Sym (3) x Xq furnished with the norm || • \\x '■= \ • | + || • ||x - 

Remark 3.1. The space of maps Xq with which we work from now on might seem 
at first glance rather ungainly. However, this space is only auxiliary in the sense that 
it does not feature in the hypotheses of Theorem 3.1 below. In our setting, it is 
essentially the 'right' space in which to capture the first correction to the leading heat 
kernel asymptotics of small initial datum solutions of the transformed equation (13). 

The main result of this section is contained in the following theorem. 

Theorem 3.1. There exists r\ > depending only on (a,b,c) G V such that for 
initial data Qq G A C H satisfying the smallness condition 

\\Qo\U<V, (14) 



12 



the Q-tensor equation (8) has a unique classically-smooth solution on M 3 x (0, oo) of 
the form 

Q(x,t) = ,4e- a *$i(x,t) + e~ at V{x,t), (15) 
where A G Sym (3) and V G Xq. 

3.1. An Auxiliary System of Integral Equations. With the above-defined Ba- 
nach spaces in mind, for any given initial data Qq G A consider the following nonlinear 
operator F : X — > X associated with Qq defined componentwise as F\ : X — > Sym (3), 
where 

r roc p 

Fi(A,V) := / Qo(y)dy+ / h(A<S> l (y,s) + V(y,s),s)dyds, (16) 

JR3 JO JR 3 

and F2 : X — > Xq, where 

F 2 (A,V;-,t) := e tA (Q -^--^f Q (y)dy] + f e^^A^s) +V(*,8),t 

rcxj r 

-*i(.,t)/ / h(A^ 1 (y,s) + V(y,s),s)dyds. (17) 

JO JR 3 

The function h is the nonlinearity in the transformed Q-tensor equation (13), and is 
given by 

h(R,t) := be- at (i? 2 - itr [R 2 ] /) - ce'^hx [R 2 ] R, 

for R G Sym (3) and t > 0. 

Suppose that there exists a fixed point (A*,V*) of the operator F in X, i.e. that 
the integral equations 

A* = F 1 (A*,V*) and V* = F 2 (A*,V*) 

hold in Sym (3) and Xq, respectively. Multiplying throughout the first equality above 
by $>i(x,t) and then adding the contribution of the second equality, one finds that the 
integral equation 

A*$i(-, t) + V*(-,t) = e iA Q + T e^ A h(A*^t + V*,s) ds 

Jo 

holds for the assumed fixed point (A*,V*). Writing R* := A*&i + V*, we note that 
the above integral equation implies that R* is classically smooth and satisfies (13) 
on M 3 x (0,oo). Following the time-dependent transformation effected by Q*(x,t) = 
e~ at R*(x, t), one immediately finds that Q* satisfies (8) on the same space-time domain. 
In turn, we may deduce from Proposition 2.1 (uniqueness of solutions) that for the 
chosen initial datum Qq G A C H, the corresponding solution Q* of the Q-tensor 
equation admits the decomposition 

Q*(x, t) = A*e- at $ 1 (x, t) + e- at V*(x, t), 

for t > as claimed. Therefore, in order to achieve this representation formula it suffices 
to show that the operator F : X — > X possesses a fixed point in X for given Qq G A. We 
achieve this for a special class of Qq (namely those which satisfy ||Qo||.4 < f] for 77 > 
sufficiently small) through an application of Banach's fixed point theorem. Establishing 
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the fact that F is locally Lipschitz on X is a little delicate, and we consider details of 
this calculation below. 

3.2. The Operator F : X — > X is Locally Lipschitz. In order to make our 
presentation more concise, we adopt the following notation. For given A,B £ Sym (3) 
and V, W £ Xq, we define the contraction modulus ko = ho (A, B;V, W) to be 

k := (\A\ + \B\ + \\V\\ Xo + ll^llxo) + (14 + \B\ + \\V\\ Xo + \\W\\ Xo ) 2 , 

and write <fi < ip whenever <p < Ckaip and C > is some absolute constant. 3 We 
demonstrate that the operator F satisfies the difference property 

\\F(A,V) - F(B,W)\\ X < \\(A,V) - (B,W)\\ X , 

for any A,B£ Sym (3) and V, W £ Xq, following which we deduce that F is indeed 
well-defined as an operator on and with range in X and thus is locally Lipschitz on X. 
Indeed, we gather these goals together in the following proposition. 

Proposition 3.2. The operator F associated with Qq £ A defined on X component- 
wise in (16) and (17) is well defined, i.e. has range in X, and also satisfies the local 
Lipschitz condition 

\\F(A,V)-F(B,W)\\ X < \\(A,V) - (B,W)\\ X , (18) 

whenever (A,V) and (B,W) belong to some bounded subset of X. 

Proof. We begin by stating a useful estimate upon which we call frequently during 
the course of our proof. 

Lemma 3.3. For A,B £ Sym (3), V, W £ X , x £ R 3 and t > 0, we have 

\h(A^x(x,t) + V(x,t),t) -h(B$i(x,t) + W(x,t),t)\ 

< (i + _^L_) 8 S e-«\t + l)- 3 \\(A,V)-(B,W)\\ x . (19) 

We now look to obtain estimates on the components of F and its differences in the 
appropriate spaces. 

3.2.1. The component F\ : X — > Sym (3). The treatment of F\ is straightforward: 
following an application of the estimate (19), we deduce that the inequality 

IF^V) - Fi(B,W)\ < dkoi^B-^WMAV) - (B,W)\\ X (20) 

holds for any A, B £ Sym (3) and V, W £ Xq and constant C\ > 0. Furthermore, one 
can quickly verify that tr[i ? i(^4, V)] = and indeed F\(A, V) < oo for any (A, V) £ X. 

We now turn our attention to the component F^, whose analysis requires consid- 
erably more care than that of F%. It proves to be prudent to cleave our estimates into 



To avoid any confusion which may arise from line to line, we interpret the string cf> < tp < u> as 
4> < Cikoip < Cik$ui. 
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one part which is in a neighbourhood of the origin in time, and one other part which is 
bounded away from zero in time. We now consider the more involved case when time 
t is bounded away from the origin. 

3.2.2. The component Fi : X — > Xq when t > 1. A judicious splitting of integrals 
faciliates a swift analysis. For (A, V), (B, W) G X, we begin by writing the difference 
as 

F 2 (A, V; x, t) - F 2 (B, V; x, t) = I(x, t) + J(x, t), 



where 



and 



/(•,*) :=jT e(*- s ) A (/i(A$i + y,s)-/i(S$i + W,s))ds, 



/•t/2 

J(;t):= J e ( *~ s)A (/i(A$i + V,s) -h(B^ + W,a))da 

poo p 

-*i(-,t) J J (KAQ^t) + V(y,r),r) - h{B^ x {y,r) +W{y,r),r)) dydr. 

We deal initially with the estimate on the term I and consider its temporal integrand, 
given by 

e —— -— (h(A^(y, s) + V(y, s), s) - h(B$i(y, s) + W(y, s), s)) dy. (21) 

[4:77(1 — s)) A i z 

for x € M 3 , t > 1 and t/2 < s < t. For any given fixed i£R 3 the inequality 

1 + -JE=) <C(1 + -#L=) (22) 

holds for all 5 > 0, t/2 < s < t and y in the region 7^i(x), where 

^i(x) := L £K 3 : < J|! 

and the constant C > depends only on (5 > 0. For this given i 6 K 3 , we consider the 
integral (21) but over the region 1Z\(x) as opposed to the whole space: 

/ $(x-y,t- s) (h(A$i +V,a) — h(B<& x + W, s)) ds 

JKi(x) 

09 \\(A,V)-(B,W))\\ X f ( \y\ y 8 - 5 

?mf Mk( 1+ N V M / Z(x-y,t-s)dy 
< \\{A,V)-(B,W))\\x ( x + Jz|_\ ~ 8 ~ S ^ (23) 



e as (s + l) 3 V v 7 ^! 
For the same fixed igl 3 , the inequality 
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holds for y in the region 1Z2{x) '■= I^ 3 \ T^i(x), whenever t/2 < s <t. Now proceeding 
in a similar manner as we did for the region IZi(x), we find that 

f <3>{x-y,t- s) {h(A$ x {y, s) + V(y, s), s) - s) + W(y, s), s)) dy 

p -\x-y\ 2 /8(t-s) 

= / 7TT+ ^e-^^ t - s Hh(A^ 1 + V,s)-h(B^ 1 + W,s))dy 

JlZ 2 (x) (4vr(t - s)) 3 / 2 

(24) / \ x \ \" 4 -| r -\x-y\ 2 /8(t-s) 

V Vt + 1/ ^ 2 (») (4vr(i - s)) 3 / 2 



e a S ( s + 1 )3 ^ ^/F+Ty Jk 2 {x) (47T(t - s)) 3 / 2 

Combining results (23) and (25) we obtain 

/ $(jc - y, t - s) s) + V(y, s), s) - s) + s), a)) dy 

£ ^ \ ^ + 4=) 4 2 ||(A^)-(i?,W))l|x. (26) 



e as (s + l) 3 V VtTT 
Finally, integrating (26) with respect to s over the interval t/2 < s < t, we find 

I(x,t)<(l + -M=) 4 2 (t + l)- 2 \\(A,V)-(B,W)\\ x , (27) 



as desired. Let us progress to the required estimate for J. 

It is at this point we utilise the heat equation 'decay improvement' Proposition 
2.2. For convenience, for our chosen (A, V), (B, W) Glwe denote by ho : M 3 x (0, oo) — > 
Sym (3) the map 

h (x, t) := h(A$i(x, t) + V(x, t),t) - h{B<5>i(x, t) + W(x, t),t). 

If we denote by M(x, t, s) the quantity 

r e -\— y\ 2 /Hts) / r \ 

M(;t,s):= — — — J -ha(y t s)dy-$(',t-8)[ h (y,s)dy) 

and also stipulate that H$ : (0, oo) — > Sym (3) be 

H o(t) := / / h (y,T)dydT, 
h Jr 3 

one may verify that the identity 



J(x, t)=J M{x, t, s) ds + H (0) t) - $i(x, t)) - $ (x, -) Ho ( 
/■V 2 / 6vr Id 2 \ 

+1 (r^-j<^0 •>*<•>* (28) 

s y / 

J4(a;,t):= 

holds for all x£l 3 and i > 1. It is straightforward to show that 

Ji(x,t)< (l + -Jl=) 4 2 ( t + ir 2 || (Ay) _ (jB , W) || X) (29) 
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for i = 3, 4 simply by using familiar properties of the heat kernel, inequality (19) and 
the fact that (t - s)- 1 < 4(i + for t > 1 and s in the interval < s < t/2. The 
treatment of the cases i = 1,2 require a little more care. An application of the result 
of Proposition 2.2 to the term J\ yields 

t/2 

M(x,t,s)ds 



o 



4+| 



(9) /"*/2 /■ | v | / 1 + \y\/\/8(t - s) , 



/0 JR3 (t - s) 2 + \ X \/y/8(t ~ s) 

Now, noticing that for < s < t/2 and t > 1 we have the string of inequalities 
(t - s)" 1 / 2 < V2V 1 / 2 < C(t + l)" 1 / 2 < C(a + ir 1/2 , we deduce from the above that 

rt/2 



/ M(x,t,, 
Jo 



-4-i 



"t> / 2 ^ I /-yi I \ fXS 



o 



x || (A, V) - (B,W))\\ X 



* 2 



< ^i + ^=j (t + l)^||(A,VO-(B,W))l|jr, (30) 

since 6 > gives us the required integrability of (1 + | • |)~ 3-5 / 2 in dimension three. It 
remains to verify such an estimate holds for the term J2, and for this we require the 
following simple lemma. 

Lemma 3.4. For x G M 3 and t > 1, i/ie heat kernel difference <&\(x,t) — $>(x,t) 
satisfies the inequality 

2 -|x| 2 /8(t+l) 

l*!^,*) - *(z,t)| < (f + 1)5/2 • (31) 
Proof of Lemma. Consider the smooth map ip defined by 

** ( l \ z 2 ( 2+t 

if>(z,t) := e 4 V2(*+D7 - e 4 l^+D 

For any fixed t > 1, by considering the equation ip z (z,t) = one may show that 
2 1 — ^ ip(z,t) is controlled by its unique global maximum over the set [l,oo), namely 
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Upon setting z = \x\ , we readily deduce from the above 



n p -\x\ 2 /8(t+l) 



from which estimate (31) quickly follows. □ 
Piecing this together with the estimates (29) and (30), we obtain 



J(x, t) < || (A, V) - (B, W) \\ x ( 1 + -J^= ) (t + l 



, - 4 -| 



which along with (27) provides 



\F 2 (A, V; x, t) - F 2 (B, W; x, t)\ < \\(A, V) - (B, W)\\ x \l + "y== J (t + V, 

(32) 

for i£l 3 and t>l. 

3.2.3. The component F 2 : X — > Xq when < t < 1. This case follows almost 
immediately from (19) and the fact that our time interval of interest in this instance is 
compact. Thus, combining this with (32), rearranging and taking norms in Xq, we find 
that 

\\F 2 (A, V) - F 2 (B, W)\\ Xo < C 2 k (A, B; V, W)\\(A, V) - (B, W)\\ x (33) 
for some constant C 2 > 0. We have therefore verified the difference property (18). 

3.2.4. The operator F 2 : X — > Xq is well defined. All that remains to be checked 
is that the map 

f e -\x-y\ 2 /4t / e -\y\ 2 /4 r \ 

(x ' t] * L -j^w [ Qoiy) " L Qo{z) dz ) dy 

is an element of Xo. The case when t > 1 may be tackled using inequality (9) and 
Lemma 3.4. On the other hand, for the case when < t < 1, we need only concern 
ourselves with the behaviour of the heat term e* A Qo near the origin in time. We obtain 
the required estimates by writing 

{e tA Q ){x)= [ $(x-y,t)Q (y)dy+ [ $(x - y,t)Q (y) dy 

JKi(x) JTl2(x) 

and applying inequalities (22) and (24) in the manner previously outlined over the 
regions TZi(x) and lZ 2 (x), respectively. 

Finally, setting (B,W) = in (33) completes the proof that F : X — > X is both 
well-defined and locally Lipschitz on X. □ 



Let 25 := B(0, eq) C X denote the closed ball of radius £o > in X, where £o is yet 
to be fixed. To close the proof of Theorem 3.1 by means of an application of Banach's 
fixed point theorem, it remains to show that F : OS — > 25 is a strictly contractive 
operator for eo chosen sufficiently small. To this end, fix eo to be the positive root of 
the quadratic equation associated with the contraction constraint 

max{C 1 ,C 2 }(2e + 4e 2 ) = ^ 
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and taking into account the definition of the contraction modulus ko, we deduce from 
(20) and (33) that 

\\F(A,V)-F(B,W)\\ X < ±\\(A,V) - (B,W)\\ X , 
for all (A, V),(B, W) G 53. Given that we also have the estimate 



\\F{A,V)\\ X < sup 1 + 



(x,t) 



4+i 



+ 



<2c 



-|-| 2 /4 



(4tt) 3 / 2 7 R 3 



+ 



Qo(y)dy 



+ -\\(A,V)\\ X , 



for any (A, V) G OS, it is clear we may find 77 > to ensure that the quantity 



sup ( 1 + 

(as,*) 



4+i 



,*A 



Qo 



,-M 2 /4 



Qo(z)dz (x) 



+ 



Qo(y)dy 



(47t) 3 / 2 7r3 

is less than or equal to eo/2 for all Qo £ A satisfying ||<2o||.A < ?7- Thus, for all such Qq 
satisfying this 'smallness' condition, the associated nonlinear operators F are strictly 
contractive with range in 03. By the Banach fixed point theorem, there exists a unique 
(relabeled) fixed point (A, V) G X and so from the discussion in section 3.1 follows the 
proof of Theorem 3.1, namely 

Q{x, t) = Ae-^Q^x, t) + e~ at V{x, t), 

for some A G Sym (3) and V G Xq satisfying \\(A, V)\\x < £o- 

Remark 3.2. We henceforth denote by ^4° the open subset of all initial data Qq G A 
for which the above decomposition result holds, namely 



A c 



#GL°°(IR 3 ) : esssup(l + |a;|) 8+5 | J R(x)| < rj 



4. The Scaling Regime L(t) = t 1 / 2 



(34) 



Before we discuss the behaviour of the correlation function c„ (r, t) as t 00, 
we comment on the behaviour of the correlation function concentrated on individual 
solutions. Utilising the decomposition formula Q(x,t) = Ae~ at <&\(x,t) + e~ at V(x,t) 
obtained in Theorem 3.1, we investigate how the quantity 



c(r,t) 



tr [Q(x +r,t)Q(x,t)] dx 



tr [Q(x,t) 2 ] dx 



behaves for large time. Firstly, noting that we have the equality 

-M 2 /s(t+i) 



$i(x + r,t)$i(x,t) dx 



1 



v / 8(4vr(t + l)) 3 / 2 ' 
and also that any V £ Xq has the norm decay property 

\\V(;t)h<C\\V\\ Xo (t+l)-V\ 
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one may then verify that 
V8e 2a *(4vr(t + 1)) 3 / 2 



tr [Q{x + r, t)Q(x, t)] dx < e -^ 2/8 ( t+1) + u(t), (35) 



trL4 2 ] 

for r£l 3 provided A 7^ 0, where the function to is bounded and continuous on [l,oo) 
and decays at least as quickly at t -1 / 2 as t — > 00. Setting r = in (35) above and 
taking the appropriate quotient, one quickly finds that 

c ( r t ) _ e ~M 2 /8t < 1 ( e -M 2 /8(t+l) _ e -\r\ 2 /8t\ _ U(t) c -\r\ 2 /8t , 

v ' ; - l + u(t) y > l + u(t) l + w(t)' 

from which we deduce 



c(r, t) — e 



= O (t~ 1/2 ) as t — ► 00. (36) 

L°°(R 3 ;dr) 



Let us reiterate that the above calculation is only valid whenever the constant matrix 
A is non-zero. If A = 0, one discovers 



c(r, t) 



tr [V(x + r,t)V(x,t)} dx 



tr 



[V(x,t) 2 ] dx 



from which no scaling information on the correlation function can be gleaned. This 
is not the case when A 7^ 0, as the map (r, t) 1— > e~' r ' 2//8t is manifestly self-similar on 
M 3 x (0,oo). We now demonstrate that the set of all initial data in A° which possibly 
give rise to solutions for which A = in the decomposition (15) is rare. 

4.1. The matrix A cannot be zero on 'large' sets of initial data in ^4°. 

It is our aim to show that the scaling behaviour of solutions is generic amongst all 
those evolving from initial data satisfying the smallness condition (14). Our notion of 
'generic' here is that the set of all such initial data in A giving rise to solutions with 
constant matrix A = constitute a closed set containing no open ball. 

Firstly, let us note that as each Qq G A° has a corresponding unique solution of 
the shape 

Q(x, t) = Ae'^^ix, t) + e- at V(x, t), (37) 
and each matrix A satisfies the explicit identity (c.f. formula (16)) 

A= Qo(y)dy+ / h{Q(y,T),T)dydT, 

we can ask about the nature of the resulting map between A° and Sym (3) which yields 
A from Qq. It is well known that real-analytic maps between open connected sets of 
Banach spaces cannot be constant on open subsets unless they are identically constant 
on the whole set. Thus, if we demonstrate the map Qq h- > A is real- analytic on A° and 
is indeed not identically the zero map, it follows that the set of all such initial data 
which produce the result A = in the decomposition (37) is rare: more precisely, it is 
closed in ^4° and contains no open ball. 
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In this direction, consider the map Fq ■ X x A — > X defined by 

T (Z(-,t),z) := e tA z + f e^ A h(Z(-, s),s) ds - Z(-,t), 

Jo 

for Z G X, z G A and t > 0, where X is the Banach space of maps 
X := |[/gC(]R 3 x [0,oo);Sym (3)) : sup ( 1 + -J^L^ ' (t + l) 3/2 \U(x, t)\ < oo 

[ (x,t)\ Vt + lJ 

endowed with the natural weighted supremum norm || • We also recall that the space 
A is endowed with the norm ||i?||.A = ess sup^. eR3 (l + \x\) s+5 \R(x)\. In what follows, 
we wield results from the theory of analytic Banach space-valued maps. We refer the 
reader to Zeidler [22] for the basic definitions and results. As a helpful first step, we 
establish the following proposition. 

Proposition 4.1. The map Fq : X x A — > X is real- analytic. 

Proof. By considering the difference 

M(R,p) + (z,0]-Fo[(R,p)], 

one may check that the first Frechet derivative dJ-$ of J-q at (R, p) S X x A satisfies 
(dT [(R,p)])(E^)=e tA ^ + bj\^ A e- as (i?3 + ER - ^tr [RE] l) ds 

-cj\^ A e- 2as (2tr[RE]R + tr[R 2 ]E) ds-E(-,t), (38) 

for all (E,^) G X x A. Owing to the fact that the nonlinearity h in the operator J-$ is 
of cubic order, the higher Frechet derivatives satisfy 

(d k J r [(R, p)]) ((Ei,£l), = in X, 

for (R,p), (Ei,£l)..., (Efc,^) G X x A and /c > 4. Furthermore, for constants C, K and 
r depending on the point (R,p), the estimates 

\\d k F [(Z,z)}\\ < — ^, whenever || (Z, z) - p)|U Xyl < r 

hold for < k < 4, where norms are taken in the appropriate space of multilinear 
operators. By triviality of the higher Frechet derivatives, the map J-q is real-analytic 
on X x A. □ 

This result aids us in establishing the useful fact that the solution map J- : Qq h-> Q 
from ^4° to X is itself a real-analytic map. 

Theorem 4.2. The solution map J- : .4° —> X is real- analytic. 

Proof. In what follows, we call (R, p) G X x A a solution pair if and only if R G 
is the unique solution of (13) corresponding to the initial datum p G A. We denote 
the corresponding solution map which yields R from p by T . We make the important 
observation that for any point (R, p) G X x A, 

To[(R, p)} =0 if and only if {R-,p) is a solution pair. (39) 
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Suppose (R, p) chosen from the open subset X x A° C X x A is & solution pair. We 
claim that the partial derivative dzJ~o[{R, p)] £ £{X) is a homeomorphism. 
As the partial Frechet derivative dzJ~o satisfies 

(d z T [(R, p)])(E) = bj* e ^ A e~ as (rS + ER- |tr [RE] l) ds 

- c j\ ( - t -^ A e~ 2as (2tr[i?H]i? + tr [R 2 ] E) ds-E(-,t). 

for all E £ X, one may then verify that the inequality 

\\d z T [(R,p)]-J\\ L{x) <1 

holds by 'smallness' in the || • ||^-norm of those solutions of (13) starting from A° initial 
data. Furthermore, the map J : X — > X given by JE := — 3 clearly has unit norm in 
C{X). It then follows from the theory of Neumann series that dzJ^oKR, P)\ e L(X) is 
a homeomorphism. Thus, by the analytic implicit function theorem, one may infer the 
existence of an open neighbourhood V C A° of p, an open neighbourhood U C X x A° 
of the point (R, p) and a real-analytic map <p : V — > X with the property that 

{(Hp),p) ■ P eV} = To 1 (o)nu. 

By observation (39), it is clear that the set J r J~ 1 (0) is simply that of all solution pairs 
{J-p, p) such that p £ A. Therefore, the map <fi : V — > X is identically equal to the 
restriction F\y of the solution operator to the set V C ^4°. 

As the point (R, p) was arbitrarily chosen, it then follows that the solution operator 
T : *A° — > X is real-analytic. □ 

By virtue of this theorem, it is clear that the map p *— > (R, p) sending an initial 
datum p £ A° to its solution pair is real-analytic on ^4°. By a routine calculation, one 
may also verify that the map 

r- r-OQ r- 

(Z,z)\-+ z(y)dy+ / h(Z(y,T),r) dydr 

JR 3 JO JR 3 

from X x A° to Sym (3) is also real analytic. With these facts in mind, one may view 
the right-hand side of the equality 

A= Qo(y)dy+ / h{Q(y,T),T)dydT 

JR. 3 JO JR 3 

as a composition of analytic maps taking X x A to Sym (3). As compositions of real 
analytic maps are themselves real analytic, we deduce that the map Qo h-> A is real 
analytic on A°. 

We need only now show that this map is non-zero on „4°. To this end, we provide 
a non-empty set of initial data contained in ^4° for which the matrix A in the decom- 
position (15) is non-zero. Together with the fact that the coefficient matrix A depends 
in a continuous manner on initial data, we have the existence of an open set of initial 
data for which 
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Proposition 4.3. For a > 0, let X a : R 3 — > R denote the map 

X «W ;= - ( i + H)^ - (40) 

If Qa,a '■= diag(A a , A a , — 2A a ), i/ten i/iere exists ao > suc/i i/iai A ^ in the decom- 
position (15) /or a// < a < ao- 

Proof. By uniqueness of solutions of (8), it follows that solutions remain of the form 
Q(x,t) = diag(A(x, i), A(x, t), — 2A(x, i)) for i > 0, where A(x,i) satisfies the nonlinear 
heat equation 

<9A 

— = AA-aA-6A 2 -6cA 3 

on M 3 x (0, oo). By Theorem 3.1, we are ensured of an e > such that 

\X(x,t)\ < Ce on M 3 x (0,oo), 

for some constant C > depending only on (a, 6, c) E P. 

Now, setting (J,(x,t) := e a 'A(a;,i), for ao > chosen small enough we have 

-be~ at (fi(x,t)) 2 - ce~ 2at (fi(x,t)f <0 on M 3 x(0,oo). 

Furthermore, since fi satisfies the equation 

fx(-,t) = e tA X - j* e^ A (b e~ as s)) 2 + 6c e~ 2as (/x(-, s)) 3 ) (is, 

by integrating across over space, we deduce 

/j,(x,t)dx< / Xo(x, t) dx < 0. 



Thus, the L 1 (M 3 )-norm of the solution cannot decay to as t — > oo, and so for the set 
of initial data 

|diag(A , A , -2A ) : A (x) = - |"|^ 8+ ^ for < a < a | C A°, 
the corresponding matrix A in the solution decomposition (15) cannot be zero. □ 

From this we conclude that the map Qq h-> A is not identically the zero map on 
A° . It follows that the set of all those Qq in A° which yield A = in (37) is closed and 
contains no open ball. 

Remark 4.1. We subsequently denote the set of all points in ^4°, for which we do 
not have information on the coefficient matrix A by B. 

4.2. Asymptotic Behaviour of c w (r, t) as t — > oo. As discussed in the intro- 
duction, we express the notion of averaging over initial conditions in a rigorous manner 
by evaluating the correlation function with respect to a suitable time-dependent family 
of Borel probability measures. We construct statistical solutions of (8) in the following 
proposition, whose proof is a modification of a construction contained in FoiAS, Man- 
ley, Rosa and Temam [8]. We quickly recall the definition a Dirac delta measure on 
the space L 2 (M 3 ). 
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Definition 4.1. For any given Q G L 2 (M 3 ), the associated map 5q : 7 ? (L 2 (]R 3 )) 
{0, 1} defined by 

f l if QeA, 

:= 

[ otherwise, 
is the Dirac delta measure concentrated at Q G L 2 ' 



Proposition 4.4. Suppose < 7 < 00. For any given .BoreZ probability measure 
~p, supported on the set 

[Q G L 2 (M 3 ) : ||Q|| 2 < 7 }ntf, 
t/iere exists a one-parameter family of Borel probability measures {/it}t>o satisfying 

*W = ^E^jSs^qJE), (41) 

/or all measurable subsets E Q K and t > 0. For eac/i k > 1, £/ie coefficients iD^j > 
respect the sum J2j$k,j = 1) S(t)Qk,j denotes the action of the semigroup of 
Proposition 2.1 on initial data Q k j. 

Proof. We establish the spaces which we use to construct the family of measures 
{l^t}t>o- We endow the set 

K := {QgL 2 (M 3 ) : ||Q|| 2 < 7} 



with the weak topology inherited from L 2 (M 3 ), with respect to which it is a compact 
topological space. Furthermore, this space is metrizable with metric d w , given by 



d w (Q,R):=jr±; [ Q: X k - [ R '■ \ 

where {x^'kLi is a countable dense subset of the set K C L 2 (E 3 ). We write Mo (A) to 
denote the set of all Borel probability measures carried by subsets of A, and furnish 
Mo(A') with the subspace weak-* topology inherited from C(K)', whence Mo(A') is 
itself a convex, compact Hausdorff topological space by the Banach-Alaoglu theorem. 
We make the observation from the estimate 



1 d 
2~dt 



which may be derived from (8), that there exists T = T(Qq,j) > such that the 
solution Q(-,t) satisfies 

Q(-,t)e {r G L 2 (1R 3 ) : ||i?|| 2 < 7 }, 

for all < t < T whenever Q(-, 0) = Qq G A n H. Furthermore, we write o~t to denote 
the set of solution trajectories 

a T ■= {Q G C([Q,T],K) : Q(x,t) solves (8) and Q(; 0) G K n H} C C([0, T]; A). 
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Using suitable properties of solutions of the Q-tensor equation (8), by the Arzela-Ascoli 
theorem ax is itself a compact topological space with respect to the topology induced 
from C{[0,T};K), where 

doo(Qi,Q2) ■= v^x i d u) (Q 1 (t),Q2(t)). 

Similarly, the space Mo(ctt) is also a convex, compact Hausdorff topological space. 

Let us now begin our construction. For a given ~p e Mq(K) supported on K n H, 
by the Krein-Milman theorem we know there exists a sequence of families of Dirac delta 
measures {&g k }j=i f° r & = 1,2,3,... satisfying 

f Nh r 

J ip dp, = lim ^f?jkj J (fdSq^ for all ip € C(K). 

Using this approximating sequence of measures, we define a new sequence of probability 
measures {yu fe }^2 =1 C Mo(ctt) by 

: = 2^^fej5g(..Q .), 

i=i 

where Q{-',Qk,j) denotes the solution trajectory starting from the initial datum Qfcj. 
This defines a sequence of points in Mo (or), so by compactness there exists a measure 
H E Mo(ot) t° which a (relabeled) subsequence of {p k }'^L 1 converges in the weak-* 
topology. We now focus on this limiting measure \x. 
For each fixed s 6 [0, T], the map 



(f^ ip(Q(s))dp 
j {j j' 

is well defined, positive and linear on C(K), so by the Riesz-Kakutani theorem there 
exists a Borel probability measure p s (depending on the choice of s 6 [0, T]) satisfying 



tp(Q(s))dp(Q) = / <p(R)diis(R) 

ax J K 

for all ip £ C(K). By direct computation, one may show for the sequence of approxi- 
mants {p k }^Li that the equality 



p(Q(s))dp k (Q) = ( ip(R)d^ k s (R) 
Jk 



holds for all ip S C(K), where 



Using the fact that fi k — fj, weakly-star in Mo(cjt), we find 

p(R)dm= [ tp(Q(t))MQ)= lim / f{Q(t))dfi k = lim / <p(R)dtf(R). 

J ax fe— >oo fc— >oo 



From this we deduce the result 

/x t (£) = lim rfiE), (42) 

fc— »oo 
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for all measurable subsets E C K. One completes the proof of the theorem by noting 
that lim^o ^t(E) = fJ-o(E) and that the transport of the measure fio may be extended 
globally in time as the semigroup {S(t)} t >o is defined globally in time. □ 

Assembling all that has come before, we now approach the proof of our main result, 
which follows rather swiftly from previous results. 

Theorem 4.5. For any given 5 > 0, there exists r] > depending only on 5 and 
the parameters (a, b, c) £ T> such that for any Borel probability measure fio supported in 
the open dense set 



A°\B= {Re L 00 ^) : esssup(l + |x|) 8+5 | J R(x)| <i]\\B, 
the associated correlation function (6) exhibits asymptotic self-similar behaviour: 



c w M)-e « 



0(t 



-l/2> 



as t — > oo. 



L°°(R 3 ,dr) 

Proof. By the previous result (41), we know that 

N k 



3=1 



as k — > oo for any measurable subset E C A° \ B. Since we have that 



C W) ( r > = , lim 



0*j / / tr [Q(x + r)Q{x)} 



dxd6 S(t)Q k iQ) 



x) I dxdS sm (Q) 



the result then follows by employing the same calculation as was used to achieve (36). 

□ 



5. The Scaling Regime L(t) = t: An Observation 

We now present evidence to support the existence of another scaling regime for the 
correlation function c w , namely L{t) = t, when is concentrated on initial profiles 
whose L 2 (M 3 )-norm is not restricted in magnitude in the manner of Theorem 3.1. 

Suppose that we offer initial data Qo £ H of the form Qq = diag(Ao, Ao, — 2Ao) for 
equation (8). By uniqueness (Proposition 2.1) we know the solution Q remains of 
diagonal form Q = diag(A, A, — 2A), where A satisfies the scalar nonlinear heat equation 

B\ 

— = AX-aX-bX 2 -6cX 3 . (43) 
dt v ; 

For any R > we denote by Aq,_r : M 3 — > R the map 




if \x\ < R, 
otherwise, 
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where A* is the global minimiser of the bulk potential associated with (43), namely 

A* = min + - A 3 + —A 4 ) . 

asm V 2 3 2 / 

In order to comment on the asymptotic behaviour of the correlation function concen- 
trated on such solutions, we require the following result (which we state without proof) 
that is based on the observations of Aronson and Weinberger [1] and Polacik 
[18] that equation (43) supports pulse-like solutions. 

Proposition 5.1. There exist R > and c > such that solutions of (43) 
subject to initial data Ao,h for R > Rq satisfy 



(44) 



where B(Q,ct) C M 3 denotes the ball of radius ct and centre 0. Furthermore, for any 
to > there exists $ = $(to) such that the solution A satisfies the bounds 

< X(x,t) < c e - ff (M- B +*) on M 3 x[t ,oo), (45) 

for some positive constants C and a. 

Utilising the results (44) and (45) above, one may verify that the correlation func- 



tion 



c Sqo (r,t) 



X(x + r, t)X(x, t) dx 
X(x, t)X(x, t) dx 



concentrated on solutions with initial data Qq = diag(Ao,_R, Aq,h, — 2Ao,_r) for R > Rq 
satisfies 

/ I T* I 



lim 

t— >oo 



C S Qo ( r ' *) - P 



t 



0, 



L°°( 



\dr) 



where P is a cubic polynomial of the shape 

P(z 



~3 (4c + z) {2c -zf 



13c 

Thus, for this one-parameter family of initial data, whose i?-norm may be made ar- 
bitrarily large, we obtain both a different universal scaling function (a polynomial, as 
opposed to the exponential map) and a different form for the coarsening length scale 
L{t). 

6. Closing Remarks 

In this article, we have been able to obtain asymptotic information for c^ (-,t) 
in L°°(M 3 ) for (jlq probing a rather restricted region of phase space H. Ideally, one 
would wish to be able to decompose phase space H = Uj^jHj in such a manner that 
given the knowledge supp/xo C Hi C H, one could find non-trivial : M 3 — > R and 
Li : (0, oo) — > (0, oo) such that 



lim 

t— >oo 



(r,t)-ri 



Li(t) 



0. 



L°°(R 3 ,dr) 
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It is in this sense we hope to find a phase portrait for c w defined on H from our earlier 
comments contained in the introduction. 

Furthermore, our main result Theorem 3.1 is a small initial data result which one 
would expect to hold for a wide class of reaction-diffusion equations, and in particular 
scalar equations. It is not clear whether or not the high-dimensional target space of the 
reaction-diffusion system (8) supports behaviour of c M distinct from that of correlation 
functions defined on scalar equations. 

Finally, it would be of some interest to investigate how the influence of flow (gov- 
erned by, say, a Navier-Stokes type system) might affect the form of the scaling function 
r or length-scale L. These are problems we hope to tackle in future work. 
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